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Abstract. We describe how to design a large class of always on spin— 1 interactions between 
polar molecules trapped in an optical lattice. The spin degrees of freedom correspond to the 
hyperfine levels of a ro-vibrational ground state molecule. Interactions are induced using a 
microwave field to mix ground states in one hyperfine manifold with the spin entangled dipole- 
dipole coupled excited states. Using multiple fields anistropic models in one, two, or three 
dimensions, can be built with tunable spatial range. An illustrative example in one dimension 
is the generalized Haldane model, which at a specific parameter has a gapped valence bond 
solid ground state. The interaction strengths are large compared to decoherence rates and 
should allow for probing the rich phase structure of strongly correlated systems, including 
dimerized and gapped phases. 
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1. Introduction 

Spin lattices are regular arrays of quantum mechanical spins with interactions involving small 
sets or neighborhoods of particles. For a lattice of spin— j particles with pairwise interactions, 
a generic pairwise coupled spin Hamiltonian is written H spm = Lj.kLa.p^j^ap^'k where the 
collective spatial indices j,k range over a D dimensional lattice and the set of Hermitian 

j2_ . 

operators {A. a }„ =0 forms a representation of the algebra u(2 j+ 1). The interaction strength 
i k 

X a p can depend both on the magnitude and the direction of the difference vector r = xj — xt 
between two lattice sites. Such models capture the essential physics of more complicated 
systems and their relevance ranges from commercial materials fabrication to fundamental 
problems in physics. Examples of the latter include studies of superconductivity and the 
origin of elementary particles and gauge fields [1]. 

Because the Hilbert space dimension grows exponentially with the number of spins, 
predicting the kinematical and dynamical properties of spin systems is challenging and has 
been the focus of decades of research [2]. Some special classes of spin lattice models can be 
solved in closed form. In those cases all properties of the associated quantum state and its 
dynamics are known. In other cases, particularly one dimensional systems [3] and some finite 
correlated higher dimensional systems [ ], the ground states have an efficient approximate 
classical description which is suitable for estimating some relevant information about the 
quantum state such as correlation functions with bounded support. In one dimensional 
quantum spin systems the dynamical evolution of a quantum state can also be efficiently 
described. By efficient, it is meant that the propagator on an n spin system can be computed 
with a number of resources that grows as a polynomial in n provided the propagation time is 
no longer than polylogarithmic in n [5]. Despite the notable progress in classical simulations, 
important properties of some spin lattices such as the behavior near a quantum phase transition 
which requires a long propagation time, or a function on the entire spectrum such as the 
free energy is required but an efficient classical algorithm is wanting. In addition, there are 
some highly correlated ground states of spin lattice Hamiltonians that serve as resources for 
quantum information processing [6]. For these reasons, it would be advantageous to be able 
to build a physical quantum mechanical simulation of the Hamiltonian. 

One approach, the stroboscopic quantum simulator, works by mapping one quantum 
system to another using a sequence of local and coupling interactions [7, 8]. Typically such 
simulations require a large overhead for spectroscopic approximation of one Hamiltonian, 
the target Hamiltonian, with another native to the simulator. Making such a procedure fault 
tolerant necessitates an even larger overhead [9]. Another approach is to build the target 
Hamiltonian directly using "always on" interactions. This works by preparing a physical 
system including a set of particles endowed with a state space isomorphic to the target spins, 
and turning on a classical control field to build the target Hamiltonian. The comparative 
advantage is a substantial reduction in control resources and time of the simulation. However, 
given that it is an analog simulation, it is not amenable to error correction and hence is not fault 
tolerant. Nevertheless, if the implemented interaction is close to the target interaction then 
many important properties of the system can be probed. Examples of atom optical proposals 
of this kind include using tunneling dynamics of neutral atoms trapped in optical lattice to 
simulate spin exchange Hamiltonians [10] and lattice gauge theories [11]. 

Prior work has proposed simulating integer spin models using neutral atoms trapped in 
optical lattices [12], [13]. There the spin is encoded in electronic ground hyperfine levels 
and the effective interaction results from tunneling induced state dependent collisions with 
nearest neighbors. A limitation of this mechanism is that the effective interaction strength U 
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is perturbative in the ratio of the tunneling rate t to the collisional interaction V, scaling like 
U ~ (t /V) 2 V. The ratio t/V must be small to maintain the state in the projected subspace 
of one particle per lattice well and imposes a fundamental limitation to the effective spin 
interaction strength. The collisional interaction itself is typically < 5kHz limited to the 
available phase space density of two atoms occupying the motional ground state of a single 
lattice well. It may be possible to build larger interactions using Feshbach resonances induced 
optically [14] or via shaped potentials [15], though it is not clear if such techniques will work 
for always on interactions. 

Recently, we proposed an always on model to generate a large class of anisotropic 
spin— 1/2 lattice models [16] using polar molecules trapped in an optical lattice. Rapid 
progress has been made in the preparation and coherent control of cold dipolar molecules 
[17] with an eye for electromagnetic trapping [18]. The spectroscopy of the molecules is 
much richer than for single atoms allowing the simultaneous trapping on optical transition 
frequencies and coherent control at microwave frequencies with weak decoherence due to 
spontaneous emission. These features make these systems experimentally relevant candidates 
for quantum simulators. In our construction the spin is encoded in the valence electron of a 2 E 
polar molecule with zero nuclear spin prepared in the ro-vibrational ground state. Spin-spin 
interactions can be induced by applying an external microwave field tuned near the transition 
between the ground and first rotational states to mix in dipole-dipole coupled states. These 
rotationally excited states are superpositions of properly symmetrized spin entangled states, 
and by judicious choice of field strength, frequency and polarization, the coupling realizes 
spin entangling interactions in the ground states. 

Here we extend that work to build integer spin models. This is a notable advance for two 
reasons. First, integer spin lattice models have qualitatively different behavior from their half- 
integer spin counterparts. In particular, antiferromagnetic Heisenberg interactions between 
one dimensional integer spin lattices exhibit a gap to excited states which is absent in the 
half-integer case. This is the celebrated Haldane gap [19]. Second, nature provides us with 
the state space isomorphic to interger spin by encoding in ground electronic hyperfine levels 
of molecules with half integer nuclear spin, of which there are many species. Engineering 
spin— j > 1/2 spin interactions is more challenging because the state space that that must be 
coherently controlled is significantly more complex. A consequence of this structure is that 
special care must be taken to negate unwanted local interactions in the designed Hamiltonians. 
Despite this, it is notable that closed form expressions for the effective interactions can be 
obtained in the asymptotic (large intermolecular separation) limit. Furthermore, our analysis 
can be used as a guide for building more complex interactions with higher spin F hyperfine 
states. 

The paper is organized as follows. In section 2 we introduce the state space of a single 
polar molecule with hyperfine structure in its low lying ro-vibrational states. Depending on the 
total nuclear spin of the molecule, half spin or integer spin can be encoded into the hyperfine 
levels. We focus on encoding in F = 1 hyperfine states and derive the asymptotic dipole 
dipole coupling between a pair of 2 E molecules, one in the first excited rotational state one 
in the ground rotational state. The effective spin interactions in the ground states are induced 
via microwave coupling to excited states and are derived in section 3. Because the interaction 
is additive, it is possible to build design a variety of spin spin interactions using several fields 
tuned near excited state potentials. We demonstrate an explicit construction of the general 
bilinear biquadratic isotropic Hamiltonian in one dimension. In section 4 we describe how one 
might experimentally probe the many body ground state of a designed Hamiltonian. Finally, 
we conclude with a summary and speculate about possible extensions of our construction. 
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Figure 1. Setup depicting a two dimensional square lattice of polar molecules. The lattice 
could correspond to an optical lattice constructed from two sets of orthogonal standing waves 
in the x and y directions and tight confinement in the y direction. The lattice is shown 
illuminated by a spatially homogeneous linearly polarized microwave field E(r). The spin- 
rotation and hyperfine couplings in each molecule are indicated. 



2. The state space: spectroscopy of dipolar molecules near ro-vibrational ground states 

The molecules we consider here are alkaline earth mono-halides with electronic ground states 
having 2 E symmetry. These molecules, like alkali atoms, have a single valence electron 
allowing for optical trapping, e.g. in a far off resonant optical lattice. Unlike the case 
with atoms, they also possess a finite dipole moment in their electronic ground states. When 
prepared in a rotational eigenstate the spatially averaged dipole moment is zero, but it can be 
made non zero by mixing rotational states with opposite parity. Because the rotational spacing 
as anharmonic it is possible to spectroscopically resolve the transition between ground and 
first excited rotational states and induce a finite dipole moment. By including nuclear spin 
on one of the constituent atoms the hyperfine interaction defines a good quantum number F 
which is the total electronic spin in the ground states. For weak dipoles induced in states with 
F = 1 the dipole-dipole interactions between pairs of molecules we thereby obtain effective 
spin— 1 Hamiltonians in the ground states. It is a straightforward matter to adapt the following 
analysis to derive effective spin models with larger integer or half integer spins. 

2.1. Single body interactions 

The Hamiltonian describing the internal and external dynamics of a pair of molecules trapped 
in wells of an optical lattice is H = H m +H ex . The Hamiltonian describing the external, or 
motional, degrees of freedom is H ex = Ya=\ Pf / (2m) + Vj(x; — X,-), where P, is the momentum 
of molecule i with mass m, and the potential generated by the optical lattice V;(x — xj) 
describes an external confinement of molecule ; about a local minimum x, with ID rms width 
zo- We assume isotropic traps that are approximately harmonic near the trap minimum with a 
vibrational spacing /ico osc . Furthermore, we assume that the molecules can be prepared in the 
motional ground state of each local potential using dissipative electromagnetic pumping [17]. 
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It is convenient to define the quantization axis z along the axis connecting the two molecules, 
%2 — Xi = Azz with Az corresponding to a multiple of the lattice spacing. We assume lattice 
trapping potentials that are strong relative to the energy scales of the dipole dipole coupled 
potentials. This means that the mechanical forces induced by coupling to the potentials do 
not significantly perturb the motional state of the ground states. The optical lattice potential 
seen by the first excited rotational state of the molecule is essentially identical to that seen by 
the ground rotational state. This is in marked contrast to the situation of atoms coupled by 
optically induced dipole dipole interactions where the excited states typically are antitrapped 
by the lattice potential. 

The Hamiltonian for the internal degrees of freedom of two molecules is H mt = 
Hdd + E/=i ff» where H m is the single body Hamiltonian Hdd is the dipole dipole coupling, 
described in detail below. For molecules without nuclear spins, H m = BN 2 +yN ■ S with N the 
dimensionless orbital angular momentum of the nuclei, and S the dimensionless electronic 
spin (assumed to be S — 1/2 in the following). Here B denotes the rotational constant 
and y is the spin-rotation coupling constant. The coupled basis of a single molecule i 
corresponding to the eigenbasis of is {\Ni,Si,Ji;Mji)} where J; =Ni+S; with eigenvalues 
E{N = 0,1/2, 1/2) = 0,E(l, 1/2, 1/2) = 2B-y, and £(1, 1/2,3/2) = 2B + y/2. Including 
effects due to nuclear spin I, the single molecule Hamiltonian is (assuming for simplicity that 
only one of the two constituent atoms has a nuclear spin) [20, 21]: 

Hi = BNj + yNj-Sj + blj-Si + cIp) 
+eQq[3lf -I(I+1)}/4I(2I-1). 

The additional terms describe respectively: the Fermi contact, dipolar spin-spin, and 
electric quadrupole couplings. The nuclear spin-rotation coupling is typically negligible 
for the 2 E molecules of interest here. A suitable basis set for the local interactions is: 
{\Nj,S,Jj,I,Fj,Mpj)} where is diagonal in the quantum numbers S,I,F,Mf but can 
couple states with different quantum numbers N and J. Expressions for these matrix element 
are given in Appendix A. We focus on the molecules where the rotational energy dominates 
over all other internal energies in H m . Therefore, we neglect couplings between different 
rotational states but include couplings (J = N + 1/2) «-> (J = N — 1 /2). Typically, the off 
diagonal matrix element is small compared to diagonal couplings meaning the electron spin 
is locked to the rotation. The spectroscopy for an examplar molecule with / = 3/2 is given in 
figure 2 

We begin our derivation of the dipole-dipole potentials, and effective Hamiltonians to 
follow, with a diagonalization of the single body Hamiltonian H m for molecules with nuclear 
spin 7 = 3/2. The rotor piece of the single body Hamiltonian 1 describes anharmonically 
spaced rotational states, so we restrict attention to the spectroscopically resolvable subspace 
N = 0, 1, The ground states with quantum numbers (N = Q,J = 1 /2,F = 1) have an energy 

E p = -S(b + c/3)/4 (2) 

A rotationally excited molecule in F — has good quantum numbers (N — l,J= 1 /2,F = 0) 
and the energy is 

E {0) = 2B+(2j-5b-c-eQq)/4. (3) 

Excited states with F = 1,2 are linear combinations of states with J =1/2 and J = 3/2. For 
a given angular momentum projection Mf, the eigenstates of H m are 

\N=1,F,±,M F ) = cos(^)|l,l/2,l/2,3/2,F,M F ) 
+ sin(^)|l,l/2,3/2,3/2,F,M F ), 
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Figure 2. Level-structure of a single 2 Z polar molecule with nuclear spin / = 3/2 within 
the lowest two rotational manifolds N = 0, 1. The spectroscopic parameters for the example 
used in this work, 40 Ca 35 Cl, are B = 4563.746MHz, y = 42.208MHz, b = 19.30134MHz, 
c = 12.4554MHz, eqQ = 1 .00284MHz, and d = 4.265 Debye [22], The bare rotational energy 
is E N =BN(N+ I). 



where the mixing angles are 



where 



!_) _ ,„„-! ( 4v/5(10fc+5c-3e q e) 
907+806- Uc+3eqQ 



<j>U,-J = tan 
(j)( 2 -) = tan 

with =(j)( R -)+7:. . The eigenenergies are 



= 2fi-(l2y-5/? + 5c)/12 
v { z ] = 2B + (30y-55b-llc-3eqQ)/60 



v 



(l0b + 5c-3eqQ)/6V5 



(5) 



(^) = <y<P+v1p)/2±J{v<P - vrV/4 + vf 2 , (6) 



v ( ! 2) = 2B-y-(b-c)/4 
vf ] = 2B+y/2-{5b + c-3eqQ)/20 
vf ] = b+(5c + eqQ)/\0. 

Excited states with F = 3 are not dipole coupled to F = 1 ground states and we do not include 
them here. 



2.2. Dipole-dipole interactions 

The near field dipole-dipole interaction between two molecules separated by r = |xi — X2I is 
H M = — ( £ (-l)*z4 9 D 2 - 9 -3Dl D20 + A.c.J. (7) 

r 9=-l 

The dipole operator coupling the ground and first rotational states of molecule i is D ; ' = 
Y} q =-\ \N = = 0,0|e*, and d is the dimensionful dipole moment. The important 
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Figure 3. Dipole-dipole coupled excited rotation states for the subspace of one shared 
rotational quantum N\+Ni = \ plotted as a function of intermolecular separation r. (a) Excited 
states for 40 Ca 35 Cl with asymptotic quantum numbers (N = 0,7 = 1/2, F;N = 1,F,±). Red 
potential curves are dipole coupled to the (0, 1/2, 1) ground states, (b) Magnification showing 
sets of potential curves used in the text for building spin interactions. This scale shows the 
relevant range of pairwise separation at optical lattice spacings. 



consideration here is that since the H^d acts to couple rotational degrees of freedom of two 
molecules, these interactions will be manifest at hyperfine energy scales due to the coupling of 
spin to rotation. Since b, c are smaller than y we have a hope to obtain tunable spin interactions 
at intermolecular separations obtainable with an optical lattice. 

The explicit dipole-dipole interaction has been solved before for the sharing of one 
optical photon between two alikali atoms with hyperfine structure [23]. After including 
nuclear spin, Hdd loses some of its symmetry with respect to operations over spin and spatial 
degrees of freedom, however it still respects conservation of total projection quantum number 
along the intermolecular axis: i.e. Mf tot — M F \ + M F 2- Additionally, [Hdd,Si.2l = where 
S i,2 is the exchange operator acting on the two molecules with eigenvalue o = ±1. In the 
subspace of one rotational quanta shared between the molecules, a good basis for H mt is 

|0, 1/2, l/2,I,Fi,M Fi ; 1, l/2,Ji,I,F!,M FH ;a) = 
■j. (\0,l/2,l/2,I,F h M Fi ) ® 11,1/2,/;,/,^^) 

+o|l, l/2 t f i ,I,Fl,M FH ) ® |0,l/2, \/2,l,F h M Fi )) . 

Denoting the set of variable quantum numbers 

S = {Fj,F[ ,Fj, Fj , Jj , Jj , M F i , M F /j , M F j , M F ij } , 

the dipole-dipole interaction in the subspace of one shared quantum of rotation is then 

\0, \/2^/2J,Fj,M Fj -,\,l/2,J'I,F'M F , r ,o) (8) 
(0, 1/2, l/2,I,Fi,M Fi ; l,l/2,j(,I,Fj,M Fli ;a\, 
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where 

A(S,I) = 6dM Fi +M Fri ,M Fj +M F ,j J [Fj , J'j , Fj , J'j] 

/_ i \F i +F j + l+2I+l+j' j +j' i +2M F j-2M F , i j J'j 1 I/ 2 \ 

1 ' \ Fi I Fj J 

/ 1 1 1/411/21/1 1 01 
\ 1/2 1/2 ^ j | Fj I Fi j\ 1/2 1/2 J> / 

Ej=-i (- 1 ) q (Fj, - M F'j\ l,q\F t , -M Fi ) 
(Fj, —Mpji 1 , -q\F!, —Mpn) - 3 (F>, -M Flj ; l,0\F u -M Fi ) 

(Fj-M Fj ;l,0\F>-M F , i ) . 

The intermolecular potentials are computed by diagonalizing H mt in blocks with good 
quantum numbers M Ftat and a. In figure 3 we plot the excited state potentials for 40 Ca 35 Cl 
with 1 = 3/2. At the intermolecular separation r « 200nm we find appreciable mixing of 
asymptotic hyperfine states. For trapping near this separation or smaller, the mixing of levels 
gives the possibility of tuning to different spin weighted states. There is also the possibility 
that the microwave fields could be tuned such that more distant molecular pairs have small 
coupling due to off resonant cancellation effects from states with energies above and below 
the field frequency. The asymptotic energy eigenvalues and degeneracies of the H^d in the 
subspace of one ground rotational molecule with F = 1 and one rotationally excited molecule 
are given in Appendix B. 



3. Effective pairwise spin-one Hamiltonians 

By applying a microwave field tuned near resonant to the dipole-dipole coupled excited state 
potentials, it is possible to engineer effective spin-spin interactions in the ground spin states 
of a pair of molecules. In the rotating wave approximation, the field interaction with the two 
molecules at positions Xi,X2 is 

2 

Hmt = - Y,( mi) i ■ efe' (kf ■ x '~° v0 /2 + h.c). (9) 

Here the electric field of amplitude Eq is characterized by the Rabi frequency |£1| = d\Eo\/h, 
polarization e F = oc_e_i + ao<?o + OC+ei, (<?o = z), and frequency 03 F . For molecules spaced 
by optical wavelengths, the dipoles are excited in phase and we can set k F ■ x, = 0. For 
convenience we choose the intermolecular axis along £. A lattice of n spins with spatial extent 
( will similarly oscillate in phase provided £\k F | <C 1, which suggests that superradiance could 
occur. In that state the spontaneous emission rate from the excited rotational state is amplified 
by a factor of n. To set the scale for these rates, note that blackbody scattering rates at room 
temperature for this system are on the order of 10~ 3 Hz compared to optical scattering rates 
induced by the confining optical lattice which are near the rate 1 Hz [18]. Smaller blackbody 
scattering rates are possible at a lower environmental temperature. Throughout, we operate 
in the low saturation limit for H m f, such that the dominate source of decoherence is from the 
optical trapping field and superradiance does not impose a strong limitation on the lattice size. 

When the saturation to the excited states is small, the effective Hamiltonian acting on the 
ground states is obtained in second order perturbation theory as 

^W=II r^-( £ „ W - £gr ) \*m> (10) 
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where are ground states with N\ = N2 = and energy E gr and {|\|/ n (r))} are 

excited eigenstates of H m with N\ + N2 = 1 and with excitation energies {E n (r) }. The reduced 
interaction in the subspace of the spin degrees of freedom is then obtained by tracing over the 
relative spatial coordinates of the two molecules: // sp i n = (H e ff(r)) Is i. In practice it is valid to 
model the interaction in the point dipole limit such that ff sp in = // e ff (Az) provided the ratio of 
the spatial extent of each trapped molecule to the lattice spacing is small (zo/Az < 0.1) (see 
Ref. [16] for details). 

In the following derivation, we ignore off resonant dipole-dipole couplings to the ground 
states with F = 2. For this approximation to be valid, we demand that effective interaction 
//spin be much smaller than the ground rotational state hyperfine splitting, or ||// S pin|| *C 
2(b + c/3), where we adopt the operator norm ||0|| = sup^j \\0 — Tr[0]ld/d||/|||\|/)||, for 
d the dimension of the Hilbert space that the operators act on (d = 9 for pairwise operators). 
If multiple microwave fields are used we further demand that Raman processes that might 
couple F = 1 to F = 2 rotational ground states are far off resonant. These terms are negligible 
for the interactions considered here. 

The effective spin interaction in the ground states as obtained from 10 is: 

He ff M = ^I£f ±) MAf(e,) ) (11) 

k=0j,± 

where the saturation amplitudes are 



s {0) (r) = m 
1 (ha F -(Ef>{r)-E st )y 

s (h±)/ r \ = h\n\(cos(^/2)-,m(¥ U±) /2)~/5) 2 

j 36(fito f -(£ J (1 ' ±) (r)-£ gr )) 

The polarization dependent coupling matrices Aj expressed in the separable basis {IgM/r) ® 
)} are given in Appendix B. For a pair of molecules separated along a rotated axis 
z' = D^(pi,p2,p3)eo, where D 7 is the spin— j Wigner rotation, the effective Hamiltonian is 
obtained using the frame rotated coupling matrices: A(e/r) — > D 1 ® D l A(D lr ef)D l1i CSi-D 1 ^. 
When multiple microwave fields are used at different frequencies and possibly different 
polarizations, the interaction is additive and complex spin interactions can be built as shown 
below. 

Many lattice models involving physical spin— j particles are described by Hamiltonians 
spanned by elements (or powers thereof) of an su(2) subalgebra of the full u(2 7 + 1) operator 
algebra. The elements of the spin— 1 irrep of su(2) are: 






In table 1 we tabulate some spin models that can be built using multiple microwave 
configurations. The quantization axis z is defined to be the intermolecular axis, assuming point 
like particles. We assume that each field is tuned near resonant to one dipole-dipole coupled 
potential such that only saturation to that state need be considered. In reality there will be 
residual off resonant couplings which will produce a Hamiltonian that deviates from the single 
state coupling result. In building many of these spin patterns, such as for the two following 
examples, some sets of potentials must be coupled to by two independent fields. This can be 
accomplished using fields tuned to both the upper and lower asymptotic manifolds. 
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Table 1. Some spin patterns that result from 1 1 . The field polarization is given with respect 
to the intermolecular axis z and the frequency of each field is chosen to be near resonant with 
one excited state potential. The interaction is // e ff( r ) — h\0\s{r)A/A with A the two body spin 
pattern operator. The saturation weights are chosen to describes the case where all constituent 
fields have the same Rabi frequency. The same spin patterns can be realized with different 
field strengths by adjusting the detunings. 
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3.0.1. Example 1 By combining three sets of spin operators in table 1 we obtain the general 
bilinear, biquadratic, isotropic interaction between a pair (j,k) of spins: 

H{* = f/(cos9S 7 - -S A . + sin0(S 7 - -S A .) 2 ), (12) 

for U > 0. This construction requires 10 microwave fields to allow for tunable parameters 
U and . We demonstrate below using numerical optimization that the model can be 
approximated with acceptable accuracy using 4 fields. 

3.0. 2. Example 2 The last spin pattern in table 1 is created using 13 fields all polarized 
along the same direction. While this is a rather complicated construction it demonstrate the 
possibility of anisotropic interactions in a plane. For a pair of spins trapped in the x — y 
plane with relative coordinate r = r(cosyy + sinyii), and with all fields polarized along £, the 
Hamlltonian is: 

H spin = U(S\S z 2 -3S]S^ n ), (13) 
where S] = S x cosy+S y i siny and U = fi|Q|(,s(r)} rel /4. 

3.1. A many-body result: the Haldane model 

We wish to extend the two body effective Hamiltonian derived above to a many body system. 
In doing so we again operate in the limit of low saturation to excited dipole-dipole coupled 
states and we treat the interaction in the ground rotational states as arising from the addition 
of the effective interaction H sp i n between molecular pairs [29]. 

A well studied spin— 1 model in one dimension is generalized Haldane model which 
involves the nearest neighbor interaction He on a chain of N spins: 

tfe=I> e /J+1 (14) 

7=1 

Ground states of this spin model have a rich phase diagram (see e.g. [24, 25]). For 
— Jt/4 < 9 < Jt/4, the system is in the Haldane phase characterized by exponentially decaying 
two body spin correlation functions and a finite energy gap to excited states that persists in 
thermodynamic limit. In that limit, the ground state is unique but for finite chains, there can 
be a degeneracy. At the point tan0 = 1 /3, each summand in 12 corresponds to a projector 
onto the spin 2 decomposition ofl<g>l=O0l02 and the 4 degenerate ground states of 
the finite linear chain are gapped valence bond solid (VBS) states [26]. VBS states have 
a practical use serving as a teleportation channel using local measurements only [27]. For 
— 3;t/4 < 9 < — 7t/4 the ground states are two fold degenerate dimerized singlets with finite 
gap to excited states. At the point = — 3jt/4, it has been conjectured [28] that there 
is a continuous phase transition to a ferromagnetic phase which for the parameter range 
7t/2 < < 57t/4 has a ground state with ferromagnetic order and gapless excitations. Finally, 
for 7t/4<0<7t/2 the systems is in a trimerized phase. 

3.2. Realization 

Consider a one dimensional lattice of trapped molecules illuminated by spatially 
homogeneous microwave coupling fields. In practice the system could consist of molecules 
trapped in a three dimensional lattice with one dimension having a much smaller lattice 
spacing relative to other two such that the dominate pairwise interactions are along that 
dimension. There is no a priori reason to expect that next nearest neighbor and longer 
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range interactions can be ignored in ID. However, by applying multiple microwave fields, 
it is possible to substantially negate the effect of longer range interactions. This occurs 
for two reasons: first, the interaction in the ground states falls off not like d 2 /P but like 
s(r) 2 (/zA — d 2 /r 3 )), so that for longer range pairs, the interaction quickly saturates to a fixed 
value; second, using additional fields turned above or below resonance, the unwanted longer 
range interaction can be negated by the additivity of H e g(r). This mechanism also allows 
for the reduction in unwanted single body interactions. We demonstrate this explicitly using 
numerical optimization to obtain a set of Hamiltonians close to Hq (see figure 4). The error of 
the implementations is measured by the operator distance of the implemented nearest neighbor 
Hamiltonian H"" to the closest isotropic Hamiltonian Hq. The field parameters for three 
specific implementations are given in table 2. The quality of the implementations do not 
change substantially using other molecular species, such as the F = 1 hyperfine levels of CaF 
with nuclear spin / = 1/2. 

Notice that the high quality implementations tend to be clustered in the first and third 
quadrants of the plane in figure 4. The reason for this bias is not completely understood, 
however, it could be partially attributed to the way we have chosen to parameterize Hq 
i.e. using a decomposition over the non orthogonal operator basis Si ■ S2 and (Si • S2) 2 . 
By comparison, we could construct an orthonomal operator basis 1.2 such that H = 
[/'(cosO'tfi +sin8'<9 2 ) where Oi =Si -S 2 /2\/3, 2 = (Si ■ S 2 + 2(Si ■ S 2 ) 2 - 81 9 /3)/2\/5, 
{/'cos 8' = Tr[(9jH] and {/'sin 8' = TrfO^//]. The probability distribution as a function of 
8' for a random Hamiltonian with fixed strength U' is p(Q')dQ' = dQ' /2k. A fixed radius 
circle in the orthonormal basis is mapped to an ellipse in the non orthonormal basis which has 
an aspect ratio [(10— 2\/T0)/(10 + 2\/T0)] 1 / 2 w 0.474 and is rotated counter-clockwise in 
the plane by an angle tan _1 ((l + vl0)/3) ~ 0.30171. Hence, absent additional structure, we 
would expect to build Hamiltonians of fixed strength concentrated along an ellipse. We hold 
to the original parameterization of 12 to stay consistent with the literature. 

Important properties of spin lattice models can be obtained from static structure factors 
which are two body correlation functions. The structure factors relative to the local operators 
Oj for a chain of N particles are defined: 

S o(<l) = ^ t eW- f HOjO),), (15) 
jj'=i 

where q assumes the role of a quasimomentum for a one dimensional lattice of spins. For 
the isotropic Hamiltonian He, relevant correlators involve the spin operator O = S z and center 
operators O = e 2mS 7/3 [25]. We have computed the correlations functions for the ground 
states of three specific implementations marked in figure 4. The ground states themselves 
were computed using imaginary time evolution of the implemented Hamiltonians over a 
matrix product representation of an infinite chain of spin— 1 particles [30]. Results for the 
correlation function Ss~(q) for an infinite chain of polar molecules (N — > ■») are plotted in 
figure 5. It is notable that even though the error of the implementations are non negligible, the 
spin correlation function is qualitatively quite close to that obtained with the closest isotropic 
model. 

4. Spectroscopic measures of many body quantum phases 

An essential component of always on quantum simulations will be measurement of properties 
of the many body state that is produced. One such observable is the energy gap between 
ground and excited states. If the system is assumed prepared in the ground state of a 
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Figure 4. Phase diagram of ground states of Hq for a chain of spin— 1 particles. The 
shaded areas correspond to regions can be prepared using engineered interactions between 
polar molecules trapped in an optical lattice. The interactions H m , H"" n are the nearest 
neighbor and next nearest neighbor interactions obtained by numerical optimization over 
frequency, intensity, and polarization of four microwave fields. The optimizations seek to 
closely approximate the desired two body nearest neighbor interaction Hq (plot (a)) while 
minimizing next nearest (plot (b)) and longer range interactions. The system parameters are set 
to those of 40 Ca 35 Cl molecules with a lattice spacing Az. = 200nm. In all the implementations, 
z polarized fields were optimal. The strength of the interaction U is constrained by demanding 
weak saturation to the excited dipole-dipole coupled states. The A, O an d EH indicate specific 
realizations with errors less than 0.05. 
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Figure 5. Spin structure factor for two implementations of nearly isotropic Hamiltonians 
on an infinite chain. Solid lines correspond to the structure factors for the ground states 
of designed Hamiltonians with nearest neighbor interactions H'™, and H™. Dashed lines 
correspond to structure factors for the ground states of the closest isotropic Hamiltonian Hq 
to each implementation. The implementation H/\ has a ground state near the Lai-Sutherland 
point [3 1 ] where period three modes of spin order are predicted. These are characterized by 
peaks in the spin structure factor at q = ±2lt/3. 



Table 2. Microwave field parameters for the realization of three implementations of the 
isotropic Hamiltonian H$ in the ground F = 1 manifold of 40 Ca 35 Cl. All fields are polarized 
along z, defined by the intermolecular axis. 



field 


frequency (Of [kHz] 


amplitude H [kHz] 


spin pattern 


1 


E^-l/h- 120.619 


18.607 


"0 


2 


/h -45.560 


19.862 




3 


EV-)/h- 188.368 


56.909 




4 


E^'^/h— 112.949 


61.667 




1 


38.456 


41.923 


Ha 


2 


E^ 0) /h- 69.242 


13.064 




3 


/A +111.868 


112.379 




4 


E^/h- 107.440 


84.709 




1 


EW/A- 32.433 


17.240 




2 


EW/h- 90.400 


25.219 




3 


E^'+VA- 64.660 


19.929 




4 


£( 0) /7i- 29.087 


23.3787 





Hamiltonian H and the excited states are coupled by single spin operators S", then the gap can 
be probed with a locally coupling field as described in Ref. [16]. Another set of observables 
are the spin structure factors 15. For our system, because of the small lattice spacing on the 
order of an optical wavelength, spin correlation functions may not be directly observable but 
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can be inferred using time of flight measurements [32, 33]. This can be accomplished by 
turning off the spin-spin interaction and the trapping fields which confine the molecules and 
allowing the particles to freely expand for a time t after which they are measured by one or 
more detectors placed in a far field zone. If the detectors can resolve different spin components 
then spatial-spin correlations are converted into momentum-spin correlations which can then 
be processed to infer relevant structure factors. Spin resolving measurement are enabled 
by using a magnetic field gradient placed in the vicinity of a detection plane far from the 
source of molecules that spatially separates different hyperfine components. This allows the 
measurement of first order correlation functions for a given spin component |oc), or second 
order correlations by observing coincidence counts for measuring a particle in state |a) at 
detector 1 and a particle in state | p) at detector 2. The latter describes a Hanbury-Brown Twiss 
type measurement and has been realized in the laboratory to measure spin independent two 
body correlation functions of atoms trapped in optical lattices [34]. There the measurement 
corresponds to coincident count measurements on pixels of a CCD camera. 

For completeness, we describe in some detail here how to infer spin dependent 
correlation functions for spin— 1 particles. The system is a chain of molecules trapped in 
an optical lattice oriented along a space fixed x axis with lattice spacing X/2 where A, is 
the frequency of the trapping light. In practice we may have an ensemble of identically 
prepared one dimensional systems trapped in a three dimensional optical lattice where the 
lattice spacing in the space fixed y — z plane is much larger than X/2 so that the chains do 
not interact with each other. For simplicity, we assume the former case so that positions 
of the source particles and detectors are along the x axis only. The latter can easily be 
treated by integrating the detected signal along a column. Consider a one dimensional 
lattice with N molecules filling N consecutive lattice sites with each molecule occupying 
the ground Wannier state centered at its local lattice well minimum. When single molecules 
are released from the trap, they freely expand into a set of plane waves. The atomic field at 
the detector plane is spanned by a set of N discrete plane waves {kj = jnX/L 2 }^_ N ^ 2 with 

mode spacing Ak = kX/L 2 , where L = y/ht/M, and M is the particle mass. The position space 
annihilation operator of a molecule with hyperfine spin component a at detector position jc/ is 
ba(xi) = L> a a (Xj)e lk ' Xi . Here a a (Xj) is the annihilation operator for a molecule in state a 
in the ground Wannier state wj of the lattice well located at position Xj = jX/2. The second 
order correlation function describing coincidence counts at detector positions x\ , X2 is then 

(2) t t 

G K a >p(xi,X2) = (ba(xi)bp(x2)bp(x2)b a (xi)) 

(aUXf)al(X e ,)ap(Xt)a a (Xj)) (16) 

+e m.2(x J -x () (at(X e )al(Xj) afi (X e )a a (Xj)}}, 

where V ~ (L 2 /47tzo) 3 is the volume of the freely evolved Wannier function in the detection 
plane. In the third line we have used the fact that (aa(Xj)ap(Xk)) = 8y.i(fla(X / -)ap(X / -)) for a 
unit filled lattice, and we have defined q\ t 2 = 2n(x\ — X2)/L 2 . For detectors that are not spin 
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resolving we trace over spin to obtain 

G (2) (xi,x 2 ) = Eo,pG^(jci,JC2) 



-^£ e ^(^^)( Sy ..S, + (SyS,) 2 }) 

where r| = 1 ( — 1 ) for bosons(fermions), Sxj,x c is the swap operator between spins S j and S^. 

In order to obtain the static structure factors 1 5 for operators diagonal in the £ basis we 
need a measurement that will output expectation values e"?u(X/-JQ) (^ a (Xt)a^(Xj)a${Xi)a a {Xj)) . 
The permutation symmetric components of these correlations can be obtained as per 1 6 with 
a = p. We cannot obtain the other terms directly using measurements in the z basis. However, 
these quantities can be obtained by measuring in a different local basis at each detector. 

First define the pairwise bases = (|oc) ± |p))/V2, and \y^) = (|a) ± ;|p))/\/2. 

Consider the q dependent terms in 16 with a,p G { — 1, 1} and a 7^ p. If we measure in the 

ct B oc B 

|xj- )i and |y ± )i bases of the detector at positional then we obtain the real and imaginary 
parts of the quantities: e m ^ x i~ Xl \a]x{X^)a\{Xj)a^{X()a^^Xj)) . Similarly, measuring in the 

|x"'' 3 )2 and ly"'' 3 )! bases of the detector at position X2 we obtain the real and imaginary parts 

of the quantities: e'i l ^ x i^ x ^ (aa(X£)al j (Xj)a^(Xe)a a (Xj)) . Combining these results with 

ct B ct B 

those obtained with both detectors are measured in the \x± )i,2, \y± }i,2 and \z)i,2 bases, 
we obtain the desired quantities e lqi ?-( x j~ x i) (^(X()ap(Xj)a a (X()ap(Xj)) . Performing these 
measurements over all pairs a, p allows the computation of operators diagonal in the spin 
basis {|m}}, and specifically the structure factors S&i (q) and Sz(q)- Local measurement basis 
changes could be done in principle by using optical Raman pulses or microwave fields to 
selectively couple pairs of states in combination with magnetic field gradients. 



5. Conclusions 



Building on previous work on designing spin lattice Hamiltonians with spin— 1/2 polar 
molecules we have shown how to build spin— 1 models using the same mechanism. The 
technique leverages off hyperfine structure in 2 E polar molecules with non zero nuclear spin, 
which is the most frequent situation found in nature. Integer spin models are qualitatively 
different from half-integer spin models and accurate simulations of such models provide 
a wealth of new physics to explore, e.g. the nature of quantum phase transitions in 
the generalized Haldane model. The results obtained here demonstrate that such models 
could be built using realistic lattice spacings (~ 250nm) with interaction strengths on the 
order of ~ 10kHz. Decoherence is dominated by spontaneous emission due to optical 
excitation by the trapping lattice. For the lattice spacing and localizations considered here the 
decoherence rate is < 1 Hz [35] suggesting a quality factor of the implementation of Q > 10 4 . 
Stronger interactions are possible with smaller lattice spacings. This could be achieved by 
eliminating the optical lattice altogether and opting for interaction induced confinement with 
self assembled dipolar crystals as proposed in Ref. [36]. 

While we have focused on optimizing one dimensional spin models, a more ambitious 
study would systematically attempt to design two dimensional spin— 1 models. Such models 
offer a host of new phenomena, among them ground states exhibiting topological order and 
emergent U(l) gauge bosons [37]. We expect that an implementation using polar molecules 
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would offer new possibilities for measurement and coherent control of such systems beyond 
that obtainable using condensed matter realizations. 
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Appendix A. Single molecule spectroscopy 

The matrix elements for the one body operator H m are [20, 21]: 

(N',l/2,J',I,F,M P \H m \N, l/2,J,I,F,M P ) = 8 NJi > 
8js (BN(N + 1 ) + \ (J(J + 1) - N(N+ 1) - 3/4)) 

+ (b + c/3) (- 1 f+l/2+J+f+I+F+l y 

/' 1/2 N \ J F J' I 



J,J',I]I(I+ 1)3/2 

+cV5(-iy' +i+F+i 



1/2 / 1 J I 1 1 J 

, , , { N N' 2 
(AT',0;2,0|iV,0)| \ \ l j U 1/2 1/2 1 

y/[N',J,J',I\I(I+l) +eQq(-l ) i/2+W+i+F-N / 4 

(N'ti-amo) \ J ' N ' 1/2 1/ F J ' 1 ' 

(/,-/;2,0|/,-/) 



(A.l) 



where we adopt the notation {•••} for the 6j and9j symbols, (ji,m\;j2,m2\j3,m^) for the 
Clebsch Gordan coefficient, and [A,B,...,Z] = (2A + 1)(2B+ 1) •••(2Z+1). 

Using the rules for angular momentum recoupling, the matrix elements of the dipole 
operator are: 



(N,l/2,J,I,F,M F \Dj\N',l/2,J',I,F',M F ,) = y/[F,J,J',N] 
{-\Y'+N+V2+i+J+f+M F ,-M F ^ F _ MF . lA ^_ MF ^ 

(n,o-, i,o\ni,o){ j f , ) x ) 2 n ;}. 

Appendix B. Asymptotic dipole-dipole potentials 

Appendix B.l. Coupling to (N = 0,F = l;N = l,F = 0) asymptotic states 
In the asymptotic regime where | |//dd('")| | ^ Y> me dipole dipole interaction is 

H ** = L ^-(-^(♦Sl + I^X^I + l^ 1 )^ 1 !), CB.D 



Designing spin— I lattice models using polar molecules 



18 



where |A£) = |0, 1/2, 1/2,3/2, 1,/n; 1,1/2,3/2, 3/2,0,0; a). The 4 excited intermolecular 
potentials and corresponding degeneracies are 



Appendix B.2. Coupling to (N = 0, F = l;iV = 1, F = 1 , ±) asymptotic states 
The 8 excited intermolecular potentials and respective degeneracies are 





- E gI - 


h£ (l.±) 


-C< 1 ' ± )(V3 + l) £ / 2 /36r 3 ; 


8i = 


1 


E^\r) 


= E gr - 


h£ (l.±) 


+ C (l ' ± \V3 + l)d 2 /36r 3 ; 


82 = 


1 


EV\r) 


Ear ~ 


h£ (l.±) 


-C( 1 - ± )2d 2 / 36 '" 3 ; 53 = 3 








= Fg r - 


h£ (l.±) 


+ C^' ± hd 2 /36r 3 ; g 4 = 3 






E^\r) 


= £ gr - 


h£ (l.±) 


-C^cP/ier 3 ; g 5 =4 








= F gr - 


h£ (l.±) 


+ C^ l ^d 2 /36r 3 ; g 6 = 4 






E^\r) 


= F gr - 


h£ (l.±) 


-d 1 ^)^- 1)^/36^; 


8i = 


1 




= F gr - 


h£ (l.±) 


+ d 1 ' ± )(V3-l)rf 2 /36r 3 ; 


8a = 


1. 



where C( L± ) = (cos(^ ±s > /2) - sin((|)( 1 ' ± )/2) V5) 2 . 

Appendix B.3. Coupling to (N = 0, F = 1; N = 1,F = 2,±) asymptotic states 
The 16 excited intermolecular potentials and corresponding degeneracies are 





= F gr - 


h£ (2,±) 


~C^ ± \l + Vl)d 2 /36r 3 - 


c?l = l 


Ff-V) 


= F gr - 


h£ (2,±) 


+ C( 2 - ± )(7 + V7)rf 2 /36r 3 ; 


82 = 1 


E^\r) 


= F gr - 


h£ (2,±) 


-C^0.251898d 2 /r 3 ; 


83 = 2 


E^ir) 




h£ (2,±) 


+ C( 2 - ± )0.257898rf 2 /r 3 ; 


84 = 2 


E?*\r) 


= F gr - 


h£ (2,±) 


-C( 2 - ± )(V3 + l)£/ 2 /12r 3 ; 


85=2 


E^\r) 


= F gr - 


h£ (2,±) 


+ C( 2 '±)(V3 + l) £ / 2 /12r 3 ; 


86 = 2 


Ef±\r) 


— Epr ~ 


h£ (2,±) 


-C^d 2 /6r 3 ; g 7 = 3 




E^\r) 


— Epr ~ 


h£ (2,±) 


+ C^ ± U 2 /6r 3 ; g 8 = 3 




E^\r) 


= Fgr - 


h£ (2,±) 


-C( 2 ' ± >0.156354t/ 2 /r 3 ; 


89 = 2 


E[ 2 ±] (r) 


= Fgr - 


h£ (2,±) 


+ C( 2 ' ± '0.156354t/ 2 /r 3 ; 


810 = 2 


E^\r) 


= Fgr - 


h£ (2,±) 


-c^Hi-^d 2 ^- 


8ii = 


E$t\r) 


= Fgr - 


h£ (2,±) 


+ C^ ± \l-Vl)d 2 /36r 3 ; 


812 = 


E^^ 


= Fgr - 


h£ (2,±) 


-C^XV^-l^/Ur 3 ; 


813 = '< 


E { ^\r) 


= Fgr - 


h£ (2,±) 


+ C (2 ^{V3-l)d 2 /l2r 3 ; 


814 =< 


Ell'%) 


= Fgr - 


h£ (2,±) 


-d 2 ' ± )0.00956771,i 2 /r 3 


815 = 


E^ir) 


— Epr ~ 


h£ (2,±) 


+ C (2 ' ±) 0.0095677W 2 /r 3 


816 = 




F gr + F(°)-2c/ 2 /9r 3 ; gi = l 
F gr + F(°) + 2,i 2 /9r 3 ; g 2 = 1 

F gr + F<°)-rf 2 /9r 3 ; g 3 =2 

F gr + F<°)+rf 2 /9r 3 ; g 4 = 2. 



(B.2) 



where C' 2,± ) = 1 — sin((f> 
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Appendix C. Catalogue of coupling matrices 

The effective spin Hamiltonians in the ground states are written in terms of coupling matrices 

(F) 

Aj weighting tensor coupling due to a field tuned near resonant with the excited state 

potential with energy Ej ' (r). These matrices are grouped according to the [N — l,F) 
excited state manifold to which that the potentials asymptote. The matrices are a function of 
the field polarization e^ = oc_e_ i + ocoeo + c*+ei and are expressed in the ground state product 
basis {| -1)| -1),|-1)|0),|-1)|1),|0)|-1),|0)|0),|0),|1),|1)|-1),|1)|0),|1)|1)}, where 
the quantization axis is the intermolecular axis. 
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